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By the example of a kicked quartic oscillator we investigate the dynamics of classically chaotic 
quantum systems with few degrees of freedom affected by persistent external noise. Stability and 
reversibility of the motion are analyzed in detail in dependence on the noise level a. The critical 
level a c (t), below which the response of the system to the noise remains weak, is studied versus 
the evolution time. In the regime with the Ehrenfest time interval so short that the classical 
Lyapunov exponential decay of the Peres fidelity does not show up the time dependence of this 
critical value is proved to be power-like. We estimate also the decoherence time after which the 
motion turns into a Markovian process. 

PACS numbers: 05.45.Mt, 05.45.Pq 



Exponential sensitivity of nonlinear classical systems 
which display chaotic behavior to arbitrarily weak per- 
turbations makes impractical the treatment of such sys- 
tems as closed ones. Inevitably, the interaction with en- 
vironment crucially influences the dynamics of such sys- 
tems. In many cases this influence can be considered as a 
noise, that turns the motion into an irreversible random 
process. 

Quite opposite, the quantum dynamics of the same sys- 
tems manifests a considerable degree of stability against 
external perturbations [l[ . A quantitative analysis shows 
that the sensitivity to an instant perturbation is of 
a threshold nature: there exists a critical value £ c of 
the strength £ of the perturbation, below which the re- 
sponse of the system remains weak. This critical value 
depends on complexity of the quantum Wigncr function, 
that can be characterized, for example, by the number of 
its Fourier harmonics. 

In general, the critical strength cr c (i) decreases expo- 
nentially within the Ehrenfest interval < t < ts during 
which the Wigner function still satisfies the classical Li- 
ouville equation and therefore the role of the noisy envi- 
ronment remains decisive. This fact is very important Q 
for understanding the quantum-classical correspondence 
in the nontrivial case of classically chaotic systems. How- 
ever, in the case when the Ehrenfest interval is so short 
that the classical exponential instability has not enough 
time to show up, in-depth information on chaotic quan- 
tum dynamics against a noisy background is quite limited 
up to now. 

In this paper we present an advanced study of this 
problem. As in [2|, we use as a typical example 



the periodically kicked quartic oscillator whose one-step 
evolution is described by the Floquet operator T = 
e -W 0) b (i^j = e-^oft+M 2 )^^), wherc 

are the bosonic ([a, a'] = 1) creation-annihilation oper- 
ators and h = a)a is the excitation number operator. 
The driving force is given by g(t) = go ^2 T S(t — r). The 
classical motion of this oscillator becomes chaotic when 
the kick strength g exceeds unity. We suppose further 
that each kick is followed by an instant perturbation 
£ T hh8{t — T) with Gaussian random intensity £ T! ((£ T ) = 
0, (£ t £t') = cr 2 (5 rr '), which models a persistent external 
noise. Such a perturbation gives rise to the phase plane 
rotation by a random angle £ T at the time moment r. For 
any given realization (history) {£} = {£1,^2, —,£t} of the 
noise, the evolution is described by the unitary operator 



u({tht) = n: 



We consider below the time evolution of the initially 
pure state p(0) = |0)(0|, where |0) is the ground eigen- 
state of the Hamiltonian H^°\ The corresponding quan- 
tum Wigner function W(a*, a; 0) = f e~Rl Q l is isotropic 
in the phase plane and occupies the phase cell H/2. Few 
(one in the case of parameters chosen below) first kicks 
produce a state of a practically general form. At a run- 
ning moment of time t > 0, the excitation of the oscillator 
and the degree of anisotropy of the Wigner function are 
characterized by the probability distributions , 



MU};i) = Hp(U};*)H 



(i) 



2 



and 



W m ({£};*) =(2-5 m0 )E 
= (2 - <5 m0 ) S 
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)-m\p({Q;t)\n) 
-o w n({(,};t)w n+m ({£};t), 
(2) 

respectively (both normalized to unity). Here t) = 

U({Cht)p(0)UH{^};t). With the noise history being 
fixed, the evolution is unitary so that the state remains 
pure during the whole time of the motion. 

Coarse-grained features of quantum evolution. 
The mean values (n){^}-t and (|wi|){£}. t calculated with 
the help of the distributions (jl!2p characterize respec- 
tively the degree of excitation and the number of 8- 
harmonics, i.e. complexity Q of the quantum state de- 
veloped by the time t. Our numerical simulations showed 
that these values do not depend on the noise history at 
a given noise level a (i.e. they are self-averaging quanti- 
ties). 

As to the distributions (jl|2p themselves, our detailed 
numerical data indicate (see upper panels in FigfTJ we fix 
system parameters as luq — 1 , ft = 1 , .go = 2 throughout 
the paper; the Ehrenfcst time < 1 in this case) that 
at a given time t they undergo (as functions of n or m, 
respectively) fluctuations, rather strong in the first case 
and much weaker in the second one, around a regular 
exponential decay with identical slopes. These slopes, 
contrary to the fluctuations, are not sensitive to the noise 
histories. Such universal exponential laws represent the 
coarse-grained distributions Q 
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which entirely ignore the fluctuations. The first moments 
(n)o- ; t and (|m|)o- ; t = (n) a -t of these distributions are the 
free parameters to be used for fitting the slopes of the 
actual distributions calculated numerically and plotted 
in Fig. [TJ Being defined in such a way, they approximate 
quite well the corresponding history-insensitive mean val- 
ues numerically found with the help of the exact distri- 
butions (|3|4|) . Moreover, the analysis of our numerical 
data reveals (see the lower panels in Fig. [1]) that the 
low moments depend very weakly on the noise level a 
as well. The chosen noise does not influence appreciably 
not only the number of harmonics but also the degree of 
excitation. Simulations with truncated excitation bases 
in the Hilbert space of different dimensions, N 3> 1 up 
to N=6000, convinced us that the rate of growth of the 
quantum mean excitation is practically indistinguishable 

2 

from the classical diffusion law (n) t = % t (straight lines 
in lower panels in Fig. [JJ). Being in contrast to the case 
of the quantum kicked rotator, this fact agrees with the 



absence of localization of the eigenstates of the Floquet 
operator T [ij]. The increase of the number of angular 
harmonics (m) t = (n) t is also linear as distinct from the 
classical exponential upgrowth. Furthermore, the ampli- 
tudes of fluctuations of the exact distributions (|H2p are 
appreciably reduced when the noise level a is growing. 

A natural way of the coarse graining consists in av- 
eraging over realizations {£} of the noise Indeed, 
as it follows from the data presented in Fig. [JJ such 
an averaging leaves the slopes practically unchanged but 
suppresses the fluctuations. The latter are entirely oblit- 
erated at the strong noise limit a — » oo. Hence we define 
finally the coarse-grained distributions as 
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(a = oo;t); 
(<r = oo; t). 



(5) 



The bar indicates averaging over the noise. Now, the con- 
nection (|m|)oo ; t = (n)oo;t directly follows from the exact 
relation @ thus leading to the identity of the slopes. 
Being plotted, the coarse-grained distributions (|3I4|) are 
indistinguishable from the distributions j5]) shown by the 
red/darkgray lines in Fig. [TJ 
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Figure 1: (Color online). (a),(b) - Probability distributions 
(11 121) at the moment t = 80 with no noise (a — 0, black), weak 
noise (a — 0.001, green/gray) and in the strong noise limit, 
(red/darkgray) In the two latter cases the data have been 
averaged over 10 3 noise realizations. (c),(d) - Time evolution 
of the moments (n)^, (\m\) a (crosses, pluses and circles for 
a = 0, 0.001 and a = oo, respectively) as compared to the 
classical diffusion law (straight magenta/gray lines). 

The information entropy. As mentioned above, 
the characteristic number (|m|) CT;t of harmonics of the 
Wigner function at a given moment of time t can serve 
0, IB] as a measure of complexity of the current quantum 
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state. Another possibility is to use Q for the same pur- 
pose the information entropy. The latter, by virtue of 
the probabilistic interpretation, can be defined as 



Averaging then over £ T , expanding up to the term ~ a 2 
and summing at last over all the moments 1 < r < t we 
obtain 
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This entropy grows monotonically with time starting 
from 1(0) = 0. Generally speaking, during the Ehren- 
fest time < t < £e the growth is linear with a slope 
determined by the classical Lyapunov exponent. After- 
wards it slows down to the logarithmic behavior. 

Sensitivity of quantum evolution to the noise. 
As usual, sensitivity of the motion to the noise could 
be characterized by the overlap ("fidelity") F({£};i) = 
Tr [p(t) p({£}; t)] of the states developed by the time t 
during the evolution with or without influence of the noisy 
environment. However, contrary to the low moments of 
the distributions {TJ [2]), the quantity defined in such a 
way is not a self-averaging one and strongly depends on 
the noise history. The appropriate reduced measure is 
obtained by averaging over all possible noise realizations 



F(a;t) = F({£};t) 



\(o\mt)\o)\ ■ (7) 



The newly defined fidelity depends on the only parame- 
ter tr instead of the entire noise history {£}. This defini- 
tion brings into consideration the averaged density ma- 
trix fM(a;t) =p({£};t) @). Since 



(n'\p^(a;r)\n) 



■ la 2 (n'-7i) 2 



(n'\fpW(o-,T-l)P\n), 
(8) 

the averaging over the noise suppresses the off-diagonal 
matrix elements and cuts down the number of harmonics 
of the corresponding averaged Wigner function. Notice 
that the normalization condition Trp( av > (a; r) = 1 holds 
during the whole evolution independently of the noise 
level. 

The unitary fidelity operator that appears in the r.h.s. 
of Eq. © reads 



/({&!*) = &(t)U({Z};t) 



-jf T ft(T) 



, (9) 



r=l 



where n(r) = T^ T nT T is the Heisenberg evolution of 
the operator n. It is easy to calculate explicitly the fi- 
delity ([7]) for the two limiting cases of weak and extremely 
strong noise. 

Weak noise limit. Keeping in the last product the 
noise intensity £ r ^ only in one certain exponential 
factor we reduce the problem to that considered in 



t— 1 m— 1 



Therefore the fidelity stays close to unity while the 
noise level a remains appreciably below the critical value 

Insert in (FigJ2|) demon- 



2/e;:> 2 >o;.. 



strates moderate decrease of fidelity during the first 100 
kicks (h ~ 1). The classical Lyapunov decay does not 
show up. This quantum regime is opposed to the classi- 
cal fast fall up to the same value in the case h = 0.01. 
Strictly speaking, the mean value (m 2 )o ;T corresponds to 
the motion with no noise. However, as it has been al- 
ready mentioned, the low moments of the harmonics dis- 
tribution are practically insensitive to the noise so that 

(m 2 >o ;T « (m 2 )^ = 2{\m\)l. T = 2(n) 2 c;r « 2 (fr)'. 

This implies that a c (t) « l/yX^i(H)^ « *~ 3/2 - 
Though the critical value decreases with time faster than 
in the case of a single instant perturbation Q, the de- 
crease is power-like as before. 
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Figure 2: (Color online) . The scaling properties of the fidelity. 
Crosses, circles and squares correspond respectively to: a = 
0.004, 0.032 and 0.256. The magenta/grey and black short- 
dashed lines show the weak and strong noise limits (|10|l and 
(|14[) ; the blue/grey long-dashed line corresponds to the fit 
(|15[) . Insert: fidelity decay in the classical Lyapunov (h — 
0.01, green/gray triangles) and quantum (h = 1, black circles) 
regimes; in the both cases a = 0, 0005. Notice logarithmic 
scale along the t axis. 



Strong noise limit. In the opposite case of extremely 
strong noise a — > oo we insert t times the completeness 
condition X)^°=o l n )( rl l = 1 while calculating the matrix 
element f(£;t) = (0|/(£; t)\0), then calculate the aver- 
age value |/(£;i)| 2 and take at last into account that 
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^ooW = E„ t |(o|^ tt |n 
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Finally we arrive at 
(n 2 |^|ni)| 2 |(ni|^|0)| 



Only the transition probabilities in the absence of noise 
are present, the forward evolution being a chain of non- 
interfering successive transitions. Quite opposite, the 
backward transition contains the entire set of possible 
quantum- mechanical paths. 

The evolution of the averaged density matrix is of spe- 
cial interest. Eq. (JSj) shows that in the strong noise limit 
the density matrix remains diagonal during the whole 
time of the motion: (n'\p ( - av ' ) (oo; r)|n) = (n'|p( d )(r)|n) = 
$n'n tci (t). The evolution equation ([8]) reduces in this 
case to 
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Figure 3: (Color online). Von Neumann entropy versus the 
time. From bottom to top (black, blue, green, magenta and 
red): a = (0.125, 1,8,64,512) • 10" 3 . Circles show the infor- 
mation Shennon entropy 1(t), see eq. ((6}. 



»P(r) = ^Q„^if»(r-1), (12) 

(13) 



n'=0 



> 0. 



The matrix Q is symmetric, positively definite, and obeys 
the condition E^/=o Qnn' = (n\n) = 1. The equation 
(fT2|) describes |9| a homogeneous Markov's chain. Notice 
that the motion does not depend in this limit on the 
properties of the unperturbed Hamiltonian H^°\ 
It immediately follows now that 



1 



n=l 



2(n) 



1 



(14) 



We have used here the exponential ansatz ([3]) for both 
w - distributions as well as the practical independence 
of the mean excitation number (n) a -t of the noise. This 
formula is in a good agreement with our numerical data. 

Moderate noise: scaling property. Analytical 
consideration is not possible in the general case of the 
moderate noise level a > er c (t). Generically, the evolution 
of the averaged density matrix p( av > (er; t) is not unitary. 
This entails state mixing and suppression of the quantum 
interference, i.e. decoherence. Nevertheless, even if this 
ratio exceeds unity, the fidelity ([7]) continues, as Fig. [3] 
clearly demonstrates, to depend only on the ratio a/a c {t) 
up to the time f(d ec )(f), when the full decoherence takes 
place and the evolution becomes Markovian. 

A simple fit (compare with 0) 



F } (a-t) 



1 



t<t {dec) (a) (15) 



describes our numerical data rather well (Fig. [5]). With 
the help of this fit, the decoherence time is estimated as 

t(dec)( a ) = \J ~§%5i where D = is the classical diffusion 
coefficient. 



Loss of memory on the initial state. The state 
mixing induced by the noise leads to loss of memory 
about the initial state. Thereupon the notion of the in- 
variant (independent of the basis) von Neumann entropy 



S(<x;t) 



Tt 



^ av \a;t)lnp iav \a;t)] (16) 



becomes relevant. This entropy increases with time 
monotonically, Fig. [3l approaching the function X(t) ([6|) 



from below when t — > t 



(dec) ■ 



After the full decoherence 



takes place, the system occupies the whole phase volume 
accessible at the running degree of excitation (n) CT:t ~ 
{n)oo-,t = (|^|)oo ; t thus reaching a sort of equilibrium. 
Henceforth the phase volume expands " adiabatically" : 
the entropy <S(c;i > tr^ec)) ~ 3-(t) m ln(|m|}oo;t remains 
practically constant when (n)^^ = (|m|)oo;t 3> 1- In par- 
ticular case of the strong noise limit, (er — > 00, t^ dec ^ — > 
0), the averaged density matrix is diagonal and the von 
Neumann entropy equals 



E^o^Wln^W 



ln(n) 



2{n) 



(t » 1) 



(17) 



We have arrived at a remarkable connection <5>oo(t) = 
T{t) . Similar connection between the information and in- 
variant von Neumann entropies has been discovered also 
in the theory of random band matrices [ltl ] . 

Reversibility versus Purity. Another interest- 
ing aspect of quantum motion under the influence of 
a noisy background is the degree of reversibility. This 
degree can naturally be measured by the mean overlap 
(Peres fidelity) of the initial state p(0) and the state 
|3(0|{£}, {£'}; t) formed during forward evolution for some 
time t under the influence of a stationary noise with the 
level er and a history {£} and then backward evolution for 
the same time under the same noise with an independent 
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history {£'}: 



¥(a;t) = Tr\p(0)p(0\{Sh it'};*) ] = 

TT[pm;t)p({e};t)}= (18) 
Tr [pM (a; t) p^ (a; t)] = P(cr; t) . 

The quantity P(cr; t) is referred for as purity Q. It can 
be expressed in terms of the mean number of harmon- 
ics (|tn|}cr;t of the averaged Wigner function. The corre- 
sponding probability distribution is related to the aver- 
aged density matrix p( av \o~,t) as 



t (a;t) = (2-6 m0 )- 



(19) 

(to be confronted with eq.©). This means in particular 
that, in the exponential approximation for mixed states 



2] 



na;t) 
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2{\m\)a;t + 1 <r^oo 



2(n) 



1 



2 In) 



1 



(20) 



Contrary to (|m|) CT;t , the mean value (|m|)cr : t strongly 
depends on the noise level and vanishes rapidly when 
<j grows. Comparison with Eq. (|14[) shows that, as in 
, the degrees of reversibility and sensitivity to external 
perturbations are directly connected, F(oo;i) = F{po\€). 
Notice also the relation S^(t) « — lnP oc (<), (t > 1) that 
follows from Eqs. (|17|20|) . 

Summary. The main goal of this Letter was to inves- 
tigate in detail the the dynamics of a classically chaotic 
quantum system with few (one in our illustrative model) 
degrees of freedom affected by a persistent external noise 
under the condition that the Ehrcnfcst time interval is so 
short that the classical-like exponential instability does 
not show up. We have shown first that the noise weakly 
influences the complexity of the quantum state, which we 
characterize by the mean number (ImD^t of 0-harmonics 
of the Wigner function. This number almost does not 
depend on the noise realization (self-averaging property) 
as well as on the level of the noise. At the same time 
the noise efficiently washes off fluctuations of the corre- 
sponding probability distribution thereby displaying the 
universal regular exponential decay of the coarse-grained 
distribution that describes the features of the motion in- 
dependent of the realization of the noise. 

The Peres fidelity that specifies a quantitative mea- 
sure of sensitivity of the motion to the noise utilizes the 
density matrix p^ av ^> (cr; t) averaged over the noise. The 



sensitivity remains weak until the noise level a exceeds 
some critical value. We have proved that with the as- 
sumptions indicated above the decrease of this critical 

value is power-like, cr c (t) w l/\jY^T=i(\ m \)oo-,T x t~ 3 ^ 2 - 
A scaling behavior has been discovered: the Peres fi- 
delity depends only on the ratio a 2 /a 2 (t) in a wide inter- 
val of the noise level up to some value <7d(i) that is con- 
siderably larger than the critical value a c (t). The scal- 
ing is destroyed only under influence of even a stronger 
noise a > ad(t). The evolution becomes Markovian in 
this case. This implies that decoherence takes place 



at the time t 



( ,U: 



h 

Da 2 ■ 



The information en- 



tropy and the invariant von Neumann entropy coincide, 
S a (t) Z(t) 7 when t > tu ec -\{cr). They are identical 
in the limit a — > oo. We have also noticed that the re- 
versibility of the motion influenced by a persistent noise 
is measured by the purity of the state at the moment of 
time reversal, ¥(a;t) = P(a;t). 
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